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We investigate the role played by the three-body DDtt dynamics on the near- 
threshold resonance X(3872) charmonium state, which is assumed to be formed by 
nonperturbative DD* dynamics. It is demonstrated that, as compared to the naive 
static-pions approximation, the imaginary parts that originate from the inclusion of 
dynamical pions reduce substantially the width from the DDtt intermediate state. 
In particular, for a resonance peaked at 0.5 MeV below the D°D*° threshold, this 
contribution to the width is reduced by about a factor of 2, and the effect of the pion 
dynamics on the width grows as long as the resonance is shifted towards the D°D°ir 
threshold. Although the physical width of the X is dominated by inelastic channels, 
our finding should still be of importance for the X line shapes in the DDtt channel 
below DD* threshold. For example, in the scattering length approximation, the 
imaginary part of the scattering length includes effects of all the pion dynamics and 
does not only stem from the D* width. Meanwhile, we find that another important 
quantity for the X phenomenology, the residue at the X pole, is weakly sensitive 
to dynamical pions. In particular, we find that the binding energy dependence of 
this quantity from the full calculation is close to that found from a model with 
pointlike DD* interactions only, consistent with earlier claims. Coupled-channel 
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effects (inclusion of the charged DD* channel) turn out to have a moderate impact 
on the results. 

I. INTRODUCTION 

Over the past decade we have witnessed fascinating progress in charmonium spectroscopy, 
especially, due to the development of 5-factories, the mass region above the open-charm 
threshold became accessible for a systematic, high-statistics experimental investigation. As 



ed "X, Y, Z states" ) with unusual 
l| • Among these new charmonium- 



a result, many new and unexpected states (the so-caf 
properties were discovered — for a recent review see Ref. 
like states the X(3872) meson found by Belle Collaboration in 2003 [2] is the best-studied 
state both experimentally and theoretically. However, the X(3872) still remains enigmatic, 
and there is no consensus on the nature of this state. 

The proximity of the X to the D°D*° threshold suggests the dynamical (molecular) 
interpretation, though other options like cc or tetraquark charmonium are discussed as well. 
The X(3872) was observed both in the J/ipTT + 7T~ (J/ipp) and J / %pTi + n~ n° (J/ipu) modes 
3|, |4j, which points to isospin violation in the wave function of the X. It is readily explained 
in the molecular picture as due to the large (about 8 MeV) mass difference between the 
charged and neutral DD* thresholds: the isospin violation is enhanced due to kinematical 
reasons, as the effective phase space available in the case of the p is much larger than that 
in the case of the u 0, Q]. A molecular interpretation implies the 1 ++ quantum numbers 
for the X(3872) and, until recently, this assignment was commonly accepted and supported 



by observation of the X in the D^D* decay mode j?-^]. However, while the analysis of the 
J/ip7r + 7T~ decay mode of the X(3872) yields either 1 ++ or 2~ + quantum numbers js], the 
recent analysis of the J/-?/>7r + 7r~7r mode seems to favour the 2 _+ assignment |4(, though the 
1 ++ option is not excluded. As shown in Ref. [loj], the X cannot be a naive cc 2 h state 
and, were the 2 h quantum numbers confirmed, very exotic explanations for the X would 
have to be invoked. In the absence of such a confirmation we stick to the most conventional 
1 ++ assignment for the X. 

Threshold affinity should lead to a significant admixture of the pertinent charmed meson 
pair in the wave function of the resonance, whatever the nature of the X(3872) is, though 
it cannot per se shed any light on the origin of binding mechanisms responsible for the 
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formation of the X. A natural explanation for the X(3872) would be a cc 2 3 Pi charmonium 
state (x'ci): residing at the D°D*° threshold, but, unless the coupling of the quark state 
to the charmed mesons channel is unnaturally small, the naive bare cc spectrum should be 
distorted strongly by coupled-channel effects. Indeed, the microscopic calculations 111. \12\ 
confirm this pattern: the X(3872) pole can be generated dynamically by a strong coupling 
of the bare x'ci state to the DD* hadronic channel, with a large admixture of the DD* 



component, see also Ref. 



13| where the fine-tuning of the x'ci state to the DD* threshold 



was discussed based on the analysis of line shapes for the X(3872). 

A competing approach is a traditional one-pion exchange (OPE) one. Historically, long 
before the charmonium revolution of 2003, pion exchange between charmed mesons was 
considered as a mechanism able to bind the isosinglet DD* mesonic system and to form 
a deuteronlike state near threshold — see, for example, Refs. Q, [l^. Immediately after 
discovery of the X(3872), the OPE model was revisited [ly, Il7| . For the most recent work on 



urther implications of 
L9| divergent integrals 



the possibility for the OPE to bind the DD* system see Refs. [1 
the nearby pion threshold are discussed in Refs. 20, [2]]]. In Refs. jl8| 
are made finite through the introduction of suitable form factors, and bound states are found 
in the static approximation for the pion and neglecting the imaginary parts of the potential. 
Only the neutral D°D*° configuration was studied in Ref. 18], and the dependence of the 
binding energy on the form factor cut off parameter A was investigated. It was shown that 
the bound state in the D°D*° system, with the binding energy around 1 MeV, exists only 
for the values of A of order of 6 GeV, that is for the values much larger than admitted 
by interpretation of the form factors in terms of quark models. The charged DD* channel 
was included in Ref. [lj]], and it was argued there that even for small cut offs of order of 
1 — 2 GeV, a bound state with a binding energy of 1 MeV appears. This result is interpreted 
then as a proof that the OPE provides enough attraction to produce a bound state. Notice, 
however, that the D*D-k coupling constant employed in calculations of Ref. 19j is too large, 
and is not compatible with the data on the D*D7r decays. 

The above-mentioned calculations treated the DD* system in a deuteronlike fashion: 
pions enter there in the form of a static potential. There is, however, an important difference 
between the deuteron and the X: the D*° mass is very close to the D°ir threshold. A natural 
worry jlj] is that in the DD* system, bound by the OPE, the pion may go on shell. The 
latter calls for the proper inclusion of the three-body DDir unitarity cuts. As shown in 
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Ref. |22j, the cut effects are of paramount importance in the charmed D a Dp system if one 
of the constituents has a large width, dominated by the S-wave Dp — > D a ir decay: bound 



states found in Ref. [23] in the static approximation disappear completely from the spectrum 
if the full three-body treatment is invoked. In the case of the X(3872) the generic two-body 
Dtt interaction goes via the D* and is in a P-wave. Due to this 'P-wave penalty" one should 
not expect disastrous consequences, though cut effects could distort strongly the resonance 
shape. Indeed, inclusion of the D* finite width alone is known to produce a spectacular 



bound-state peak in the D°D°7i° mass distribution — see Refs. (24jj27|. This justifies a full 
investigation of the role played by the three-body dynamics on a near-threshold resonance, 
which is the subject of the present paper. 

One of the most important findings of our study is that the X-dynamics, and especially 
the value of the effective coupling constant X — > DD*, is, in the molecular scenario, com- 
pletely fixed by the X binding energy Eg, as long as Eb "C AM, with AM = 8.08 MeV 
being the distance to the next (charged) two-body threshold, in line with the properties of a 
true two-body state, although there are various thresholds near by. This is in line with the 



results of Ref. 



28], where the interplay of scales was studied for the case of the presence of 



5-wave interactions only. At the same time, by an explicit calculation, we have shown the 



validity of the central assumption underlying the X-EFT [2l|, |29|, namely, that pion effects 
can be treated perturbatively for most observables. At leading order in such an EFT, pions 
can be integrated out and predictions for the observables can be made based on universal 
asymptotic behaviour of the DD* wave function 301 ]. 

The most striking effect of dynamical pions is seen in their effect on the imaginary parts 
from the DDir intermediate states. Specifically, the part of the X width stemming from 
the width of the D* gets cut in half once dynamical pions are included. A similar effect is 
observed once the three-body DDtt cut is accounted in the DD* potential: the imaginary 
part of the DD* potential gets reduced by a factor of more than 2 compared to that for the 
static pion potential which leads to a further reduction of the X width. In total the width 
of the X from the DDtt intermediate states is reduced from 102 to 44 keV due to the effects 
of dynamical pions. This observation could be of relevance for the X line shapes below the 
elastic threshold. 

The paper is organised as follows. In Sees. II A and B we introduce the notations 
and derive the system of dynamical Faddeev-type equations for the DD* scattering. In 
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FIG. 1: Kinematics of the DD* scattering due to the OPE. Double lines denote Z)*'s while single 
lines denote -D's. 

Sec. II C we apply these equations to the X(3872) assuming 1 ++ assignment for this state. 
Physical implications are discussed in Sec. III. First, we compare the residues of the DD* 
scattering amplitude calculated with purely contact DD* potential (with no pions) to our 
full dynamical results. We also discuss the role of coupled- channel dynamics for the residue. 
Then, assuming that the X(3872) has a resonance state with the peak at Eb = 0.5 MeV 
below the D°D*° threshold, we calculate the D D n line shape within our fully dynamical 
treatment and compare it with various approximations. The summary of the most important 
results is given in Sec. IV. 

II. THREE-BODY FORMALISM 

A. Kinematics and main definitions 

Consider first three channels: 

|2) = DD*, \2) = DD*, \3)=DDtt, (1) 

coupled by the OPE — see Fig. [T] with m' = m and m'^ = since, for a time being, we 
ignore isospins; the isospin structure of the interaction will be considered in detail below. 
Here m and m* are the mass of the D(D) and the bare mass of the D*(D*), respectively. 
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The dynamical equations that emerge are basically a diagrammatic representation of 



those presented in Ref. 



311 ] . In the centre-of-mass frame the momenta in the two-body 



systems DD* and DD* are defined as 

Pd = P, Pd* = ~P, (2) 
Pd = Pd* = -P, (3) 

while, in the three-body DDir system, the momenta can be defined in terms of two sets of 
Jacobi variables, {p, q} and {p, q}: 

tti tti 
Pd = P, Pd = -Q : P, Pn = q ^— P, (4) 

or 

Tfl TYI 

Pd = -q ; P, Pd = P, Pw = Q - 1 — P (5) 

m + m + m n 

Jacobi variables belonging to different sets are related to each other as 

q = aq + f3p, p = -q - ap, (6) 
q = aq + (3p, p=-q- ap, (7) 



where 



m o 2 1 ( 2m + mir)™^ 

a = — , P = a -1 = — — — 

m + m w [m + m^) 



The DD*tt vertex is 

v D d*M = 9 e • q, (8) 

where e is the D* polarisation vector, q is the relative momentum in the Dir system, and g 
is the coupling constant — in general here one could also introduce a form factor, however, 
since we are interested in near-threshold phenomena, such a form factor would not change 
the general properties of the results. The constant g can be fixed from the D*° — > D°7r° 
width via 

r(D*° -> D°7i°) =r m = ^g 2 fi q (D°n ) [2/i g (£> 7r )(m * - m - m n0 )} 3/2 , (9) 
where the reduced mass is defined as 

(vv\ m xm Y ( . 

^ XY ) = ~ , rn ■ ( 10 ) 
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Then the channels ([I]) communicate via the interaction potentials 

V£(p,q;p') = gq m 5(p-p'), 
V£(p,q;p') = gq m 5(p-p'), 

and similarly for and V%?. 

The inverse two- and three-body free propagators are defined by 

2 2 2 

D 2 (p) = m + m ^ + ^-M, D 3 (p, q) = 2m + m T + i^ + ^— 



(11) 
(12) 



(13) 



with M being the mass of the system, which can be expressed, for example, in terms of the 
energy related to the DD* threshold, while the reduced masses are defined as 

mm* m{m + m^) , „ s mm„ 



m(m + m n ) , . 

Hp = — ■ , fX q = H q {D7T) 

2m + m n 



m + m* 

and, nonrelativistically, /i* = fi p . 

Finally, we define the bare Dn self-energy as 



m + m-n 



(14) 



E(p) 



9' 



q 2 d 3 q 
D 3 (p,q) 



(15) 



B. Three body equation for the t matrix: One-pion exchange potential 



The dynamical equation for the t matrix t^ n (p, p'), where i, k = 2, 2, 3 and m and n are 
the Lorentz indices, reads schematically: 



t = V- vg t, 



(16) 



where Qq is the diagonal matrix of the free Green's functions (fT3|) and the interaction po- 
tential V possesses just two nontrivial components given by Eqs. ( TTTT) and ( fi~2l) . The system 
of nine equations (TIB]) splits into three decoupled subsystems: 



(17) 
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(18) 



^22 — V32 77"^22; 
L>2 
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/ 1 

^23 — ^23 — ^23—^33 
T 

t2 3 = ^23 - V^—tga (19) 

1 1 

^33 — —^32—^23 — V32— ^23- 

Below, in this chapter, we give schematically the derivation of the three-body equations 
for the t matrix which follows from Eqs. ( TT71) - (|T9l) . To simplify notations distinguish neither 
between V23 and V23 nor between V32 and V32. We also omit all arguments. The full and 
detailed derivation can be found in the Appendix. 

In the systems (fTTI) and ({TBI we exclude the third equation and, using the definition (fl5l) . 
arrive at 

1 1 1 

^22 — — E + S— t 2 2 + ^23 77" ^32 77-^22 

I 1 1 1 ( 20 ) 

^22 = — ^2377" V32 + E77- + ^2377" ^3277-^22, 
JJ 3 LJ 2 JJ 3 L>2 

and a similar system for the components t 2 2 an d ^22! which we do not quote here. Note that 
the relations between the breakup amplitudes £23 and t§3 an d the two-body amplitudes t 2 2 
and ^22 are given in Appendix A. 

The interaction respects C parity, so it is possible to define C-even and C-odd DD* 
matrix elements. Indeed, the system (120]) can be rewritten for the combinations 



t 2 2±*22, (21) 



which satisfy then the following equations: 



A-^ti = -E T ^377-^32 ± ^377-^3277-^, (22) 

L> 2 U3 U 3 L>2 

where the inverse dressed D* propagator A(p) is introduced as 

2 

A(p)=m* +m + ^ M-E(p). (23) 

Finally, substituting 

^ = -^ + X a± X' (24) 

one arrives at the following equation for the new function a™ n (p, p'): 

a± = V±- ^±A- X a ± , (25) 
which, in the full form (see the Appendix), reads 

aT(p,p\E) = Vr(p,p) ~ f d 3 sVr P (p,s)A- 1 (s)a p ± n (s,p',E), (26) 
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where the generic OPE potential has the form: 



D 3 {p,p') m^ + m 

In what follows we confine ourselves to the C-even states only and therefore we consider 
only the amplitude a™ n (p, p') = a mn (p,p') and the corresponding potential V™ n (p, p') = 
V mn (p,p'). 

Equation (1261) was derived neglecting isospin. However, in the case of the X(3872), the 
two two-body thresholds, the neutral D°D*° + c.c. threshold and the charged one D + D*~ + 
c.c., are split by only around 8 MeV, so they are both potentially relevant, and Eq. ff26l) is 
to be modified accordingly. 

For antimesons, we stick to the convention of Ref. [ljj, that is we define 

\M) = C\M), (28) 

where C is the C-parity transformation operator, and \M) (|M)) denotes the wave function 
of the meson (antimeson). Then interpolating currents of the pions and the D and D* meson 
involved are chosen as 

tt° = —^=(uj 5 u - dj 5 d), 7r + = ir~ = uj 5 d, (29) 

V2 

D° = uj 5 c, D° = cj 5 u, D + = (J75C, D~ = cj 5 d, (30) 

D*° = u^c, D*° = -c^u, D* + = d^c, D*- = -c^d. (31) 

Since the pion is an isovector, the isospin structure of the OPE potential ( )27|) is given by 
the product T\-t%. Therefore, in order to evaluate the isospin coefficients, we are to compute 
the matrix element {i\f\ ■ f^l A:) for i, k — 0, 0, c, and c, where the latter states are defined 
as: 

|0> = D°D*°, |0) = D°D*°, |c) = D + D*-, |c) = D~ D* + . (32) 

This will define the full OPE potentials V^ n (p,p'). Notice that, since strong interactions 
respect C-parity, then the following relation holds: 

■ T 2 \k) = (i|fi • T 2 \k), i,k = 0,c. (33) 

It is easy to find then for the two nonvanishing coefficients: 

(0|ti ■ r 2 |0) = (c\t\ ■ r 2 |c) = 1 (neutral pion exchange), 
(0|ti ■ t 2 \c) = (c|ti ■ t 2 \0) = 2 (charged pion exchange). 
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It is clear therefore that nonvanishing OPE potentials are (see Fig. [I]): 

V$r(p,p') = V^ n (p,p'), with m = m' = m , m w = m^o, (35) 

V™(p,p') = V™ n (p,p') 7 with m = m' = m c , m n = m,o, (36) 

V™(p,p') = V£ n (p,p'), with m = m , m' = m c , m w = m w c, (37) 

VoT(p>pO = Vg^ n (p,p'), with m = m c , rri = m , m w = m w c, (38) 

where m , m c , m T o, and m^c are the masses of the neutral and charged D meson, and the 
pions, respectively. The factor of 2 needed for charge-exchange potentials [see Eq. (jMj) ] will 
be introduced in the scattering equations explicitly. 

The explicit form of the OPE potentials can be obtained then as the generalisation of 
the generic potential (1271) : 

VrM) = (p' + oc ik p) m {p + a' ih prF ik {p,p'), F ih (p,p') = - - (39) 



D 3 ik(p,p') 



with 



, ^0 _ , _ m c 

^oo — a oo — ; > a cc — OL — ■ . 

m T o + mo m n o + m c 

, m c , m 



a 0c = « c0 = : > "cO = « 0c 



and 



m^c + m c m^c + mo ' 



£>3oo(p,p') = 2m + m,o + -^ + i- + V - M - zO 



p 2 


p' 2 
2m 


(P + P') 2 


2m 


2m 7r o 


p 2 

2m c 


-f h • 

2m c 


(P + p'f 


2m n o 



D 3cc (p, p') = 2m c + m,o + ^— + ^— + ^ ^ 7 — M — iO, (40) 

2m c 2m x o 

D 30c (p,p') = m c + m + m n c + ^— + ^— + ^ - M - iO, 

2m 2m c 2m 7r c 

V 2 T>' 2 (v + p') 2 
D 3c0 (p,p') = m c + m + m n c + + % - M - iO. 

The 16-component OPE interaction potential V^™" (p,p') gives rise to a system of coupled 
equations for the 16 components of the DD* scattering t matrix. However, due to symmetries 
of the OPE potential (!35|) - (l38|) . many t matrix components coincide with one another. In 
particular, it is easy to demonstrate that 

«r = a-r> ar = «3r> a = o, c , (4i) 

so that only eight independent t matrix components remain and split into two groups, 
four components in each, which satisfy two disentangled subsystems of equations. In what 
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follows we shall be interested in the |0) final state, so we consider only the first subsystem 
of equations of these two, which reads 



( r d 3 s 

a™(p,p>,E) = V T(p,p') - J ^-^V™ p (p,s)a^(s,p',E) 

~ 2 J -^K P (P,s)a^(s,p',E) 

a™(p,p',E) = 2V™(p,p') - 2 J -^V™ p (p,s)a%(s,p',E) 



(42) 



V™v(p,s)a%(s,p',E). 



X( 

\ 

d 3 s 

C. 1 ++ channel 

From now on we stick to the quantum numbers inherent to the X(3872) charmonium, 
assuming the latter to be 1 ++ as per discussion in the Introduction. To perform the partial- 
wave decomposition of the scattering equation ( l4"2l) we expand the amplitude and the po- 
tential in terms of the spherical vectors Y JLM (n): 

a r(p,p'i?) = EE fl * Ll ' i2 (^' E )E(^^ n )) m ( y V/K)) n . ( 43 ) 
vrM) = EE vi MM M E^^Hn^M^or, (44) 

J L1L2 M 

where n and n' are the unit vectors for the momenta p and p', respectively, and, given 
the chosen quantum numbers, the relevant matrix elements are those with J = 1 and 
Lx,L 2 = 0,2. To simplify notations, we omit everywhere the superscript for the total 
momentum J — 1. 

Furthermore, once we are interested only in the S'-wave in final state, we need to know 
only the af k s and af k s matrix elements of the amplitude, so that we use the following de- 
composition of the amplitude: 

aT( P , P', E) = af k s (p, p\ E)Tgg + af k s (p, p\ E)T™, (45) 

where, using properties of the spherical vectors, one can find for the projectors: 

Tgs = ~~^_^mm T DS = — — y=(5 mn — 3n m n m ) (46) 
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and, consequently, for the four nonvanishing matrix elements of the potential: 

V$ {P,P') = y J F ik(p,p\x) (a ik p 2 + a'i k p' 2 + (a ik a' ik + l)pp'x) dx, 



v£ d (p,p') 



I jFifc ^' P '' ^ i^' ikP ' 2 + aikp2 ~ 2~) + ^ aika ' ik + l )PP >: ^j dx > 



3 



~~3~ J Flk ( p > p '> x "> (^teP 2 + a 'ikP 2 Q^ 2 - + i a ika'ik + l)pp'x^j dx, 
V£ D (p,J) = 2n f\ k {p,p',x) Q(a lfc p 2 + a' lk p") ^ - ±) + + - pp'x 



3 / ( 5 3 3 



where x = cos 9, with 9 being the angle between p and p' . Then, finally, we arrive at the 
system of four coupled equations: 

a s 00 s (p,p',E) = V* 8 {p,p') 

sMs -V^{p,s)a^{s,p',E)- [ 4^Voo D (P,s)ag s (s, P ',E) 



• S2 ' /,S V f ( P ,,)afo ? ( S , P ',E)-2 / ^V£?(p,s)a° s (s,j/,E) 



A c (s) uc ^' ' cu v " ' ' J A c {s 
a° s (p,p',E) = V° s (p,p') 

S ds T _ n e / >. cc / / _ x /* 5 ds 



V D s (p,s)a s 00 s (s,p',E)- / ^ttf D (P,*)^W,-E) 



A (s) uu " ' uuv ' f ' ' 7 A ( S ) 

- 2 / £^yC S (P^)4 5 (^p'^)-2 / £^ D (p,s)a? 5 (^£) 
a s J(p,p',E) = 2V 8S {p,p') (48) 

- 2 j^V 8S (p,s)a s 00 s (s,p',E)-2 J £^V c r( P ,s)a° s (s,p',E) 

s2dS -V^ip,s)a^(s,p',E)- [ 4^Vc C D (p,s)a° s (s,p>,E) 



A c (s) - ' J A c (» 

a° s (p,p',E) = 2V c f (p,j/) 

- 2 J^V° s (p,s)a s 00 s (s,p>,E)-2 J ^V° D (p, s)a° s (s,p> , E) 
s2< ' S V° s (j>, S )a s J(s,p>,E)- [ 4^V° D (p,s)a° s (s,p>,E). 



We choose to define the energy relative to the neutral two-body threshold: 

M = m * +m + E. (49) 

The full inverse propagators A and A c entering the system of equations (1481) are obtained 
as a generalisation of Eq. ff23l through the introduction of the running width r(p) which 
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incorporates both the effect of the self-energy as well as contributions from other D* 
decay channels: 

2 2 

Ao(p) = m * + m + -^ M-Vo(p), A c (p) = m c * + m c + M-V c (p), (50) 

where the reduced masses are 



The running widths -To(p) and -T c (p) take the form 



(51) 



r (p) = r(D*° -+ D 0f y) 



(52) 



57T 



2/x 9 (D u 7r u ) M - m - m - m w o 



2/i, 



3/2 



x M — m — m — m w o 



2/i, 



-i^(£>V) 



2/i q (D 7v°) I m + h m + m,o - M 



P 



2/2 



0* 



3/2 



x ( mo + h m + m T o — M 

2/i * 



+2^( J D+7T-) 



p. 



2fi q (D + n ) M - m — 77 m c - r?v 



2/i 



0* 



3/2 



x M — m — m c — m^c 



2/i, 



0* 



-2i[i q (D + n-) 



2/j, q (D + 7i ) I m + \- m c + m n c — M 

2/i * 



3/2 



x ( m + h m c + m n c — M 

2/iQ* 
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8tt 
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(53) 
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FIG. 2: The amplitude for the D°D°7r° 5-wave production rate. The cross denotes the pointlike 
source T . 
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Note, formally there should be a radiative term also for the decay of the charged D*, however, 
in this case it is negligibly small and will therefore be dropped. 



The D D°ir S"-wave rate is calculated now as (see Fig. [5]) 
dBr B 



where 



I /-'max 

<IE 2 7 r / , / ,(.2 / ,, / (.D<'7r'U;) ; V2./ pdp J p _ pdp 



xl E + E*- 



jf 
2^ 



J(p,E) 



E + E* 



P 

2fi % 



J(p,E) 



MP) 



i 

+- 

a 



(a 2 + 1){E + E.) 



J(p,E) = l 



p2 _|_ p2 
2/ip 



oo „2 



Re 



J(p,E) (J{p,E) 



A (p) V A (p) 



s 2 <is 

o A^y 



a^(s,p,E), 



P 



max,min 



'2(x q (D°Ti ) [E + E, 



2/ip 



± ap 



(54) 
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m (m + m n o) m 
£/* = m * - mo — m^o, /x p = — ■ , a = a 00 - 



2mo + m n o m n o + m 

The overall coefficient B = |J-"| 2 absorbs all the details of the short-ranged dynamics respon- 
sible for the X production. The last term in the curly brackets in Eq. (154p corresponds to 
the interference of the production amplitudes with the pion produced by the D°* and D°*. 
In agreement with earlier claims 30[, for the bound state, interference affects substantially 
the magnitude of the production rate below the D°D°* threshold. 

Finally, for numerical calculations, we use the following masses and widths: 

m = m{D°) = 1864.84 MeV, m c = m(L> ± ) = 1869.62 MeV, 

m * = m(D*°) = 2006.97 MeV, m c * = m{D* ± ) = 2010.27 MeV, (55) 

m n0 = m (vrO) = 134.98 MeV, = m( 7 r ± ) = 139.57 MeV, 

j; = r(D*° -> D°n°) = 42 keV, r(D*° ->■ D 7 ) = 21 keV. (56) 

Then, using Eq. ((9]), the coupling constant g can be extracted from the D*° — > D°ir 
width to be: 

g = 1.29 • 10~ 5 MeV" 3/2 . (57) 



D. Regularisation and renormalisation of the three-body equation 

The system (|4"8]) is to be solved numerically. Because of the P-wave D*Dir vertex [see 
Eq. (|8])], the integrals on the right-hand side diverge linearly. We separate the short- and 
long-range dynamics of the system, where the long-range interaction is due to the OPE while 
the short-range one, in addition to the short-range part of the OPE, may contain other 
contributions, for example, due to a strong coupling of the DD* to the quark-antiquark 
charmonium. This short-range dynamics is parametrised by a constant (in momentum 
space) Co (A) (below referred to as contact interaction or as counter term) which appears 
just as an extra term in the potential V£ s , that is now 

27r f 1 

Vf k S (P, P') = Co(A) + — J ^ F ik (p, p', x) (a ikP 2 + a' ikP ' 2 + (a lk a' ik + l)pp'x) dx, (58) 

while all other components of the potential are left intact. Here A is the cut off parameter 
which regularises the integrals in Eq. (T48l) (we use the simplest sharp cut off prescription, 
that is we substitute J °° ds — > J^ds). Ideally one chooses for every A the value of the 
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counter term in such a way that it neutralises the dependence of the physical observables 
near threshold on the cut off A. As we will discuss below, this does not work for some values 
of A at least for the particular regularisation used to solve the system (1481) . This may indicate 
the need of additional counter terms. However, since the residual A dependence is very mild 
in a large range of cut offs (300 MeV < A < 1700 MeV and 2500 MeV < A < 3800 MeV), 
the renormalisation with a single counter term appears to be appropriate to investigate 
the issues at hand — the complete discussion of the properties of the system P"8"]) will be 
postponed to a subsequent publication. 

III. THE EFFECT OF DYNAMICAL PIONS ON A NEAR-THRESHOLD 

RESONANCE 

We will now compare the solution of the full problem ( 14*81) (including the counter term 
discussed in the previous chapter), with the solution of the same problem in the static 
approximation. This implies using the static OPE potential and omitting the momentum 
dependence of the D*° width, cf. Sec.[ffl""B]and Eq. ([79]) for more details. In both calculations 
we require the system to possess a resonance at E = —Eb with Eb = 0.5 MeV. Note that 
in the presence of pions there is no bound state, so we define the position of the resonance 
as the position of a peak in the D D n production rate, given by Eq. (I54"l) . Therefore, 
comparing the properties of the results of the two calculations we investigate the role played 
by dynamical pions. In addition we study the role played by the charged channel in the 
formation of the X(3872) state. 

A. Investigation of the DD* contact interaction 

Before we proceed, we briefly discuss the calculations for a model with the contact DD* 
interaction only [just retaining the Cq term in Eq. (I58l) ]. The advantage of this model is that 
it is solvable analytically and thus a comparison of this simplified model with the results of 
the full calculation allows one to understand better the findings. In the next subsections the 
corresponding equations for the single-channel as well as for the coupled- channels problem 
are derived, for the explicit solution will allow us to better understand some properties of 
the full equations. 
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1. Single- channel case with a contact DD* interaction 

We start from the simplest, single- channel, case of Eq. ( 14*8]) with the only nonvanishing 
potential being V^ 5 = Co. Then Eq. ( HHl) reduces to a simple algebraic equation: 

a oo = Co — Coa 00 I , (59) 

where, in order to simplify notations, we set in this chapter a 00 = a^Q . Here 

C A s 2 ( k ( A - k\\ 

h{E)= ds^— , _ = 2/x * A + - iTr + log— — , (60) 



.,„ s 2 /(2// *) - E - iO n V 2 \ °A + fc y/ 
with A; 2 = 2/xo*E. 

The strength of the contact interaction Co is fixed by the requirement that the system 
possesses a bound state at E — —Eb (Eb > 0): 

Cq" 1 = -Io(-Eb), aoo 1 = I (E) - I {-E B ). (61) 

To proceed we remind the reader that Eq. (I59p is nothing but a toy model for the full system 
(|48p where the cut off A is to be chosen large enough, that is larger than the natural scales 
in the full dynamical problem. This implies that A ^> s/JIq^Eb at least, so that the integral 
Iq can be expanded in powers of the ratio k/A, 

i , k 2 \ /A; 4 



I « 2/io, (A + -nk - - \ + O ^ , (62) 
which corresponds to the effective-range expansion of the scattering amplitude: 

a oo = -/^Q*^ ( _a_1 ~ ik + ^ refffc2 ) ' ( 63 ) 

where the scattering length and the effective range take the form: 

11 4 

r cS = -v- (64) 



y/2^E B 1 - y/2nQ*E B r c ft/2' ' 7rA 
The amplitude a 00 possesses two poles in the complex k plane: 

ki = i^2fi *E B , k 2 = i(— a/ 2 /^*-^ J • 

W / 

In particular, for E B = 0.5 MeV and A = 500 MeV, one finds 



(65) 



a = 6.6 fm, r cfT = 0.5 fm, k x = z31.1 MeV, k 2 = i754.3 MeV. (66) 
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Finally, expanding k at the bound-state energy for \E + E B \ <C E B , 

E + Eb 



k = y/2/i 0il: E ps iy/2nQ*E B ^1 



we arrive at the expansion of the amplitude a Q0 at the pole: 



(67) 



ReS Q 00 _ gg / fi oN 

ao ° - = ~eTe~b~' m 



where the residue reads 



Res a 00 = — 1 tKo(Eb) = ~^—tz /0 F v 69 



with k (Eb) = y/2jj *EB- This residue squared is shown in Fig. [3] for A = 500 MeV as the 
red (upper) dotted line. Its mild deviation from the straight line, at least for E B < 3 MeV, 
is caused by a nonzero but small value of the effective range, so that we have, approximately, 

g 2 cS oc ^Eb~- (70) 

Summarising, one can see that we have a large scattering length and one near-threshold 
pole in the k plane. The residue squared (that is g^ s ) of the amplitude at the bound- 
state pole scales almost linearly with the bound-state energy. These results allow for a 
nice physical interpretation using the method suggested by S. Weinberg in the mid 60's 



321 ] . He showed at the example of the deuteron, that an analysis of low-energy observables 
allowed one to quantify, in terms of the probability factor Z (0 < Z < 1), the admixture of 
the bare ("elementary") state in the wave function of a near-threshold bound state. This 



approach was generalised in Ref. 33j to the case when inelastic channels are present, as well 



as to the case of an above-threshold resonance. In particular, in Ref. 33], the connection 



is established between the admixture of a bare state and the structure o 
singularities of the scattering amplitude. The latter was considered in Ref. 
resonance classification. The results of Refs. 



' near-threshold 



33 



34| as a tool for 



34j can be briefly summarised as follows. A 



state is mostly elementary, if there are two nearby poles in the scattering amplitude, which 
corresponds to a small scattering length, and large and negative effective range. In this case 
the coupling of the state to the hadronic channel is small, and the resonance line shape takes 
a Breit-Wigner form. The state is mostly composite, if the scattering length is large, and 
there is only one near-threshold pole. The effective range is of natural size (of the order 
of range of forces) and plays a role of a correction. This scenario requires a large coupling 
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FIG. 3: Residue of the D D scattering amplitude squared versus the binding energy in the D D 
system. The upper, red (lower, blue) dotted curve corresponds to the solution of the single(two)- 
channel D°D*° problem with the contact DD* interaction. Solutions of the full three-body equation 
with dynamical pions are given by the solid lines: upper, red line — for the single-channel case and 
lower, blue line — for the two-channel case. The straight dot-dashed line (black) is shown to guide 
the eye. 

of the state to the hadronic channel. In this case the wave function is dominated by its 
molecular component. The residue of the relevant pole is determined by the binding energy 
of a composite particle and thus it is model-independent. For the parameters used as an 
illustrative example above [see Eq. (1661) ] we clearly have a predominantly molecular state. 



In the case of two coupled channels, the scattering equation for the two components of 
the amplitude, for the sake of simplicity denoted in this chapter as aoo = clqq an d a co = af(f , 
takes the form: 



2. Two-channel case with contact DD* interactions 




(71) 



with the solution for aoo given by 



Ooo — 




(72) 



(l + Co/o)(l + Co/ c )-4C 2 / / c ' 
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Similarly to the single-channel case, the integrals are 

= jf " tmJ-E-n K 2 "»' ( A + -2^ - T ) + (I) ' (73) 
= f " s/^ + 'L-E-n * ( A + & - 1) + (§) • < 74 > 

where fcg = 2/i *£, /c 2 = 2// c *(£ - AM), and AM = m* + m c - m* - m = 8.08 MeV. As 
before, we require that there is a bound state at E = —E B , that is the amplitude (1721) has a 
pole at this point. This leads to the equation to determine the strength of the contact term 
Cq. Notice that now, in the two-channel case, this equation is quadratic, as opposed to the 
linear equation in the single- channel case: 

(l + C h{-E B )) (l + C I C (-£ B )) - AClh{-E B )I c {-E B ) = 0. (75) 

It is easy to verify that the latter equation always possesses two opposite-sign solutions 
for the Co- In particular, for A = 500 MeV and E B = 0.5 MeV these solutions are Co = 
13.442 x 10" 7 MeV" 2 and C = -4.425 x 10~ 7 MeV" 2 . 

Notice that the coupled-channel problem is still dominated by one relevant model- 
independent near-threshold pole which, for E B = 0.5 MeV and A = 500 MeV, is still located 
at k\ = z31.1 MeV equal to the value of k\ for the single-channel case given in Eq. (IB~6l ). 
The position of other poles is very model-dependent, however they all are in general located 
quite far away from the threshold. 

Following the same lines as in the single- channel case above, one can find the residue of 
the amplitude a 00 at E = —E B . The dependence of the residue squared on the binding 
energy for A = 500 MeV is illustrated in Fig. [3] by the lower (blue) dotted line. To a good 
approximation, the dependence of the residue on the binding energy can be written as 



const (A) k c (E b )k q (E b ) ( 8 k c (E b )kq(E b ) x 

V<>h '" iin k c (E b ) + k (E b ) { L ttA k c (E b ) + k (E b ) I - ' l,) 



where k c {E b ) = ^J2^{E B + AM), k (E b ) = \/2/j, *E b , and const(A) is a constant of order 
of unity. In the limit of small binding energies (E B AM), one has hq(E b ) k c (E b ), 
and the formula for the residue (J76l) can be simplified as 

Res a o ~ const ■ °^ ^ ( 1 ^-n (E B ) J , (77) 

TTHfr \ ttA J 

so that, similarly to the single-channel case [compare with Eq. (|69|) ]. the residue squared 

remains approximately linear with E B , up to small finite-range corrections. In principle, the 
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range of validity of this expression is rather limited since the correction ~ \[Eb from k c {Eb) 
starts to play a role quite rapidly. However, an analogous correction but with an opposite 
sign appears also from the finite-range term (~ 1/A) in the denominator. This leads to a 
partial cancellation of these terms. 

Interestingly, for Eb ~ AM, the linear behaviour of the residue squared as a function 
of Eb is still approximately preserved since again the finite-range corrections are cancelled 
to a large extent by the corrections due to k c (Eb)- As a result, for A =500 MeV, the 
behaviour of the residue squared is almost exactly linear in the whole range of Eb- This 
linearity is, of course, accidental — for larger cut offs the deviation from the straight line is 
more pronounced; however, it remains a correction to the leading linear behaviour. 

Deviations from the Weinberg prediction are expected to be of order r^/2/i *-E_e, with 
r being the range of forces. If we estimate the latter by r e ff from above, 
0.1a/ Es/MeV, we expect a deviation of the residue from the Weinberg prediction of the 
order of 30% at Eb = 7 MeV, in line with what can be read off from Fig. [3] for the single- 
channel curve. As explained above, accidental cancellations make the result for the coupled- 
channel calculation more consistent with a linear behaviour than expected naively. 

The influence of the coupled-channel dynamics on the residue can be better understood 
if one considers the limit A — > oo. In this follows from Eqs. fl69l) and (1761) . the ratio 

of the residue for the coupled-channel problem to that for the single- channel one reads 

^ Res aoo (coupled-channel) I j E~b 

Res aoo (single-channel) a->oo I y Eb + AM 

where we neglected the tiny difference between /i * and \i c * and used the fact that const (A) — > 
1 when A — > oo. It is clear from the ratio (ITS]) that the coupled-channel effect is negligible 
at small binding energies, Eb <C AM, that is, once we stay close to one threshold, its 
effect dominates over the effect of the other, remote, threshold. In the opposite limit, for 
Eb ^> AM, when the splitting between thresholds can be neglected, this ratio tends to 
one half, that is, in agreement with natural expectations, asymptotically both thresholds 
contribute equally to the residue. 
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3. Counter-term in presence of dynamical pions 

The behaviour of the counter term Co versus A in the full single-channel problem (D°D *) 
with dynamical pions is shown in the left panel of Fig. H] by the blue dotted curve. The 
figure demonstrates a clear limit- cycle behaviour of the contact term with the increase of 
A in full analogy with the NN 



351 and 3N 



problems, see also Refs. 



37 



for related 



works. Meanwhile, a plateau between the first negative and positive infinite solutions is 



about 7 GeV for the DD* problem which is much larger than that for NN 35]. 

The behaviour of the counter term Cq on the cut off A in the two-channel case, when 
dynamical pions are present in the problem, is illustrated in the right panel of Fig. HI We 
therefore solve the full problem — see Sec JII CI In contradistinction to the one-channel case 
(and to the NN case), in a coupled-channel problem there are two solutions for Co (A) that 
correspond to a resonance with the peak at E B = 0.5 MeV [see the quadratic equation 
( I75p and the discussion in the previous section]. These solutions, shown as solid red and 
dashed black curves in the right panel of Fig. HI again exhibit a sort of limit-cycle behaviour 
individually. However, it is interesting to observe that there are regions of cut offs where 
no solution for Co(A) exists. This can be understood as follows: for values of A, where the 
two solutions for Co(A) approach each other (see, for example, the region of A just above 
2 GeV), a second singularity approaches the DD* threshold. As a consequence, the two 
approaching states start to repel each other to avoid a crossing of the levels. Therefore there 
is a region of Eb where the peak of the resonance cannot be reached by any variation of 
the strength of the potential Co (A). The regions of A corresponding to the absence of the 
resonance peak at Eb = 0.5 MeV for any values of Co(A) are illustrated in the right panel 
of Fig. HI Further details on the limit cycle behaviour in many-channel problems will be 
provided in a subsequent publication. 



What should be stressed here is that as long as we choose the cut off to be relatively 
far away from the problematic region one of the states is always located far away from the 
threshold and therefore should not affect observables. As a consequence, observables are 
basically independent of whether one or the other branch of Co is chosen — see Sec lHI Bl 

A comment is in order here concerning the role played by the OPE for the binding of the 
X meson. From Fig. H] one can see that for particular values of the cut off (for example, for 
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FIG. 4: Left panel: Behaviour of the counter term Co versus A in the full single-channel problem 
(D°D *) with dynamical pions. Right panel: Behaviour of Co as a function of the cut off A in the 
full two-channel problem with dynamical pions. 

A « 900 MeV in the coupled-channel problem), there is a solution with Co = which could 
be naively taken for the bound state in the pure OPE potential with dynamical pions. This is 
not the case however since the details of the short-ranged dynamics are simply hidden in the 
particular cut off pattern for the divergent integrals taking place for such A's. Furthermore, 
in previous chapters, we showed how to arrive at the bound state with a purely contact 
potential. Finally, the bound state can be found in a combined short-range+OPE potential 
with Co varying from — oo to oo, as shown in Fig. HI Since physical observables, for example, 
the width of the X-meson (see below) are almost A-independent in a large interval of A's, 
then no model-independent statement is possible concerning the importance of the one-pion 
exchange for the formation of the X(3872), contrary to earlier claims 
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19] 



B. The effect of dynamical pions on a near-threshold resonance 

With the experience gained in the previous chapter, we are in a position to perform a 
detailed numerical analysis of the system (148]) . We compare the following three situations: 

1. The single-channel problem in the static approximation, that is we solve the system 
(gHD with V™ n (p, p') = V™ n {p, p') = V™ n (p, p') = and with the static OPE potential 
which, as before, is given by Eq. (139|) . but with a o = « o = 1 anc ^ with 

£foo ic (P> P') = 2 ^o + m,o + {P + P ' )2 - (m * + m ) - iQ (79) 
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for the inverse three-body propagator. In addition, instead of the running width (|52|) . 
we use a constant width 

T = 63 keV. (80) 

Note that in the absence of the three-body dynamics both these effects, the constant 
width and the imaginary part of the potential, represent the same effect of the decay 
D* — > Dn related to two-body unitarity. Thus, keeping only one of these effects would 
lead to inconsistent treatment. 

2. The full dynamical calculation for the single-channel problem, including the three- 
body %DD intermediate states as well as the dynamical width of the D*. For the D*° 
width we use formula f[52j) without the contribution of the charged channels. 

3. The full dynamical calculation of the two-channel problem, including three-body DDtt 
intermediate states as well as the dynamical width of the D*, as was explained in 
SecJITCl 



First, we investigate the effect of dynamical pions on the residue of the scattering am- 
plitude aoo squared — see Fig. [3] and the discussion in Sec lHI Al The upper (red) and lower 
(blue) solid curves correspond to case 2 and case 3 above, respectively. As was discussed 
before, the solid curves lie on top of the dotted curves indicating that the binding energy 
dependence of the residue is basically independent of the dynamics. The closeness of the 
results to the straight line (black dot-dashed curve in Fig. [3]) indicates the dominance of the 
single-pole scenario and thus the molecular nature of the X(3872) meson. 

In Fig. [5]we show the results for the production rate dBr/dE corresponding to all three 
different cases described above: (i) solution of the single-channel problem in the static limit— 
(green) dot-dashed line; (ii) solution of the single-channel dynamical calculation — (blue) 
dashed line; (iii) solution of the full two-channel dynamical problem — (red) solid line. In all 
cases the value of the cut off A was fixed to 500 MeV and the strength of the contact operator 
Co was adjusted such that the D°D*° scattering amplitude has a resonance state peaked 
at E = —Eb = —0.5 MeV. All curves are normalised to obey the same energy behaviour 
near the three-body D D tt° threshold, located at E = —7.15 MeV. The difference between 
them is visible only in the region of the peak and is pronounced in the different strength of 
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FIG. 5: Production rate (in logarithmic scale) for the three calculations as described in the text: (i) 
solution of the single-channel problem in the static limit — (green) dot-dashed line; (ii) solution of 
the single-channel dynamical calculation — (blue) dashed line; (hi) solution of the full two-channel 
dynamical problem — (red) solid line. All curves are normalised near the D°D°ir threshold, located 
at E = —7 MeV. The inlay shows a zoom into the peak region in linear scale. 

the peaks — see Fig. |5j If one normalises all three curves at the peak, the difference in the 
resonance region shows up only in the different widths of the resonance. To quantify the 
effect we notice that, in the energy region around the resonance peak, the DD* scattering 
amplitude aoo acquires basically a separable form: 

1 



a^(p,p',E) 



/i0*71" 



-f{E)(j>{p)4>{p') 



(81) 



with cf) being a formfactor, so that the energy-dependent amplitude f(E) is factored out from 
the integration over the phase space in Eq. (j5l|) . Parametrising f(E) by the Breit-Wigner 
(BW) shape, one finds for the rate: 

/ dBr\ const 

\dEj BW = (e + E B y + ri/4 



kes(E), 



(82) 
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where k eS (E) corresponds to the integral over the phase space in Eq. (15^|) with J(p, E) = 1. 
Except for the region very close to the threshold, k e s(E) is a very smooth function of energy 



around the peak as compared to the Breit-Wigner shape. Therefore, expression ( 1821) can 
be used to extract the width i~x- The results for i~x extracted for all three cases above are 
shown in Table [fl Examining the transition from case 2 to case 3, one can see that including 
the charge channel changes the shape of the spectrum only marginally The nearest threshold 
which affects the observables due to the coupled-channel effect is D + D°7r~ (and D~D°tt + ) 
which is about 2.5 MeV above the D°D*° threshold. Also the D + D~tt° threshold is just 
3 MeV above the D°D*° threshold. Had the mass of the charged D-meson been smaller by 
2.5 — 3 MeV and more, the new cuts would have occurred in the DD* potential increasing 
the X(3872) width and thus emphasising the role of the coupled-channel dynamics. In the 
meantime, the transition from case 1 to case 2 elucidates the role of the three-body dynamics: 
switching on three-body effects makes the resonance narrower by a factor of 2. 

In order to quantify the role of the pion dynamics in the vicinity of the threshold we 
use the scattering length approximation and study the impact of the three-body dynamics 
on the D°D*° scattering length. The molecule scenario for the X(3872) implies that, near 
threshold, the effective-range corrections are small (which is also in agreement with the 
nearly linear behaviour of the residue as a function of the binding energy in Fig. [3]), so we 
resort to the scattering length (SL) approximation for the production rate: 



A similar approach to the analysis of the DDir line shape was used in Ref. [25], where the 
width of the X originated from a constant D°* width in the absence of inelastic channels. In 
Eq. (IBB"]) we also neglect inelastic effects from remote thresholds and introduce 72 to account 
for the change of the X width due to the DDir intermediate states in the energy region 
near the DD* threshold. Stated differently, 72 accounts for the shift of the pole position 
in the complex plane due to the three-body dynamics. Parameters 71 and 72 (72 > from 
unitarity) corresponding to the full dynamical treatment can be well determined from a fit 
to the resonance structure in the line shape — see Fig. [5j 





1 



(83) 



27 



Case E B , MeV r x , keV 71, MeV 72, MeV 
Case 1 0.5 102 31.1 1.58 
Case 2 0.5 53 31.1 0.82 
Case 3 0.5 44 31.1 0.68 



TABLE I: Parameters of the distributions (|82p and (|83|) extracted from the fit to the peak for the 
line shapes depicted in Fig. [5j 

Near the peak, that is for \E + Eb\ <C Eb, the scattering length formula ( 1831) is identical 
to the Breit-Wigner one fl82l) with the following identification of the parameters: 

E B ^tpi, r x = 2 -™i. (84) 

2/^0* f^o* 

Furthermore, for a narrow resonance (i~x <C -Eb), these relations can be inverted to give: 



^v 7 ^^, 72^ ff^ - (85) 

Clearly, a reduction of the X width by a factor of 2 would necessarily suppress 72 by the 
same amount. The results for 71 and 72 for three different scenarios from above are quoted 
in Table HI It should be noted that the expression fl83|) provides a very good approximation 
for the exact line shape in the near-threshold region even relatively far away from the pole, 
for \E + Eb\ ~ Eb- 

We would like to stress also that the effect of pion dynamics depends strongly on the 
position of the resonance state relative to the relevant two- and three-body thresholds. 
Indeed, once it approaches the three-body Z) D 7r° threshold, the imaginary parts originating 
from the three-body effects vanish, and the X width tends to a constant due to the D°* — > 
D°7 transition, while the X width in the static approximation stays nearly constant — see 
left panel of Fig. |6j On the other hand, the dynamical width of the X(3872) grows as the 
resonance peak approaches the two-body D°D*° threshold. The dependence of the effective 
X width i~x [as well as that of the 72 by virtue of Eq. (185]) ] on the binding energy can be 
parametrised in the interval from 0.1 to 2 MeV in an analytic form: 



r x{E B ) = -^- 2 + -^- 2 , (86) 
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FIG. 6: The X width as a function of the binding energy Eb- Left panel: The results of three 
calculations as explained in the text; the notation of curves is as in Fig. [5j Right panel: Comparison 
of our results [red solid (dashed) curve including (omitting) the D*Dj width] with those of Refs. [21! 



30] (dot-dashed line for the LO calculation of Ref. 

n 

Ref. [21|], respectively). 



301 ] and blue band for the NLO calculation of 



with r x = 35.5 keV, ft = 2.369 MeV, F 2 = 28.7 keV, and ft = 0.364 MeV. In this form the 
effect of dynamical pions can be conveniently incorporated in the analysis of experimental 
data. 



In the right panel of Fi g. [6] the results of our full dynamical calculation are confronted 



with those of the V-EFT 
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3Cj. In leading order (LO) of V-EFT |30| pion degrees of 



freedom are integrated out, and the wave function of the X, as a bound state, is taken in a 
universal form provided by the contact DD* interaction. Then the X width is calculated as 
an integral over the three-body phase space from the product of the X wave function and 
the decay vertex D*° — > D°ir°. The result of Ref. 30] is improved in Ref. 2l| by taking 



next-to-leading (NLO) corrections. The counting scheme used in Ref. [2l[ treats the pion 
effects perturbatively in distinction from the nonperturbative inclusion of the LO contact 
operators in analogy to what was proposed for the NN system in Refs. 4o|. The results 
of these two calculations (taken from Fig. 7 of Ref. 2l| ) are shown in the right panel of 
Fig. [6] by the dot-dashed curve (LO) and by the blue band (NLO), respectively. For the 
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sake of comparison our result 1 is given in the right panel of Fig. [6] as the (red) solid and 
dashed lines, where the former is the result for the full calculation, while the latter is the 
result when the DDj channel is omitted — since this contribution to the width was not 
included in the X-EFT calculations, it is this curve that is to be compared to the results of 



Refs. 



21 



30(. We are therefore to conclude that the X-EFT reproduces the results of our 
more complete calculation very nicely. This level of agreement provides strong support of 
the power counting underlying the X-EFT. 

All the results shown in Figs. |5] and |6] were calculated with A = 500 MeV, and the value 
of Co for the full problem (case 3) was chosen to be negative, see solid (red) branch in the 
right panel of Fig. HJ Note, however, that the variation of the results with the cutoff is very 
mild for the range of A from 300 to 1700 MeV and from 2500 to 3800 MeV. For example, 
varying A in the range 300 MeV < A < 1700 MeV influences the width only at a few percent 
level comparable with the numerical noise. Also we obtain that choosing another branch of 
Co results in the uncertainty in the width of a similar (small) size. 

As the final remark let us mention that, aiming at the effects related to dynamical pions, 
we considered only the part of the X(3872) width that comes from the DDn intermediate 
state. The physical X(3872) supposedly acquires a significant part of its width from inelastic 
channels like J/iptnr and J/ipTTim. However, the effects discussed here should be of relevance 
for the line shapes of the X(3872) in the DDir channel. 

IV. SUMMARY 

We investigated the role played by dynamical pions on the structure of the X(3872) 
by solving the full three-body DDir problem. The two-body 7rD — > ttD input is fixed 
by the width of the D*, however, a short-ranged DD* — > DD* contact term needed to 
be introduced to arrive at well defined equations. Depending on the regulator used for 
solving the equations, the contact term can take values between — oo and +oo for the cut 



1 We would like to stress that in both above-mentioned calculations the X is assumed to be a genuine 
bound (not a resonance) state. In the full treatment, however, such a bound state does not exist, so that 



one cannot calculate the width directly, following the lines of Refs. 
the X width using the procedure described above in this chapter. 
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30j . Instead, one is left to extract 
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off regularisation scheme. On the other hand, the strength of the one-pion exchange is 
fixed. Thus, the relative importance of the pion exchange as compared to the short-ranged 
contributions to the DD* scattering potential appears to be strongly regulator-dependent. 
We are therefore led to conclude that no model independent statement on the importance 
of the one-pion exchange for the formation of the A(3872) is possible, contrary to earlier 



claims 



16 



19|. 



In addition, we found that the residue for X — > DD* is weakly dependent on the kind of 
pion dynamics included. Especially, the dependence of the residue on the A binding energy 
is very close for a fully dynamical calculation and for a calculation with a contact-type 
interaction only. A deviation between the coupled-channel and the single- channel treatment 
is clearly observed but with the larger effect for binding energies beyond 1 MeV. 

The most striking effect of dynamical pions is observed in their impact on the X line 
shapes: in the fully dynamical calculation the width from the DDtt intermediate states 
appears to be reduced by about a factor of 2, from 102 down to 44 keV, assuming that the 
A(3872) corresponds to a resonance state with a peak at 0.5 MeV below the DD* threshold. 
Stated differently, by using the naive static approximation for the DDk intermediate states 
one overestimates substantially their effect on the X width. Although the total width of 
the X is rather saturated by inelastic channels like J/ifmn and J/ifminr, our findings should 
be of relevance for the predicted line shape in the D°D°7r° channel, where the signal below 
the DD* threshold is controlled exactly by the DDti cuts. For example, if the line shape 
near threshold is analysed within the scattering length approximation, the imaginary part 
of the inversed scattering length, 72, which is the only parameter affected by the three-body 
dynamics, has to be taken from the fully dynamical calculation. The parameter 72 (and thus 
the role of pion dynamics) appears to be a smooth but vivid function of the binding energy, 
for which we provide a simple analytic parameterisation suitable to mimic three-body effects 
numerically in the data analysis. 

On the contrary, the effect of the coupled-channel dynamics on the X width turned 
out to be rather moderate, which can be attributed to the fact that both the real part of 
the resonance pole E% and the X width Ix are small as compared to the separation AM 
between the neutral and the charged thresholds. 

Our results for the X width and for the residue appear to be in good agreement with those 
obtained in the A-EFT approach 2l|, |30], which justifies the central assumption underlying 
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the X-EFT that pions can be treated perturbatively. In addition, an important progress 
achieved in our work due to considering the full coupled-channel dynamics with nonper- 
turbative pions is that pion range corrections are included to all orders in our calculation. 
Because of the proximity of the DDtt threshold, those are expected to give a prominent 
contribution to the range corrections even beyond leading order. We therefore expect the 
nonperturbative calculation to have a smaller uncertainty than that with the perturbative 
treatment of pions. In addition, the X pole is located very close to both the DDtt as well 
as the DD* thresholds, which all influence the pertinent integrals. Thus, although for the 
X the effects can be nicely absorbed into a counter term, this is not necessarily the case 
anymore for other, related resonances, say in the B sector 2 . In this sense our calculation 
appears as an important step for a common understanding of a larger class of quarkonium 
resonances. 
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Appendix A: Details of the derivation of the three body equation (1261) 



After excluding the third equation, the systems (117)1 and ( ITS]) , with all arguments and 
sub (super) scripts restored, read 

E(p) 



t™(p,p',E) = -Z(p)5 mn 5(p-p' 

+ g 2 f d\ 



2 / j3 wi 

(aq p + f3pp) 



D 3 (p, q)D 2 (-q - ap) 



t™(p,p',E) 



*22 VP,P,E) = -g — \- —— 

D%{p, —ap — p) D 2 (p 

3 Qm 

(aq p + (3p p ) 



■q?{p, P \E) 



+ 9 Z d'q 



D 3 {p, q)D 2 {-q - ap) 

E(p) 



t%(-q-ap,p',E), 



tW(P,P',E) 



t^(p,p',E) = -Z(p)5 mn 5(p-p') , 

2 [ j3- qm{aq p + PPp) pn/ / 

+ 9 I d q ^ - (-q -ap,p, E) 



DsiP, q)D 2 {-q- ap) 



t^(p,p',E) = -g 



2 (ap m + p' m )(ap' n +p r . 



+ 9 Z d*q 



D 3 (p, -ap- p') _ D 2 {p) 

3 _ q m (aq p + (3p p ) 



^(P) +mn{- -I j?\ 
^T*22 KPiPiE) 



t^(-q-ap,p',E). 



l D 3 (p,q)D 2 (-q-ap) 22 
It can be demonstrated then that t^(p] p', q') and tl^(p; p', q') of the form 



t™(p;p',q';E) = gq'J{p - p') 



- t2 T(p,p',E)q' n 



(Al) 



(A2) 



L>2{-q - OLp) 

f£(p; p', q', E) = g4j[p - p') - ^Ly^P, P', E)q' n 

^ 2 n (p -q' - ap>, E)(aq' n + (3p' n ). 



L>2{-q - olp 

satisfy the last system, given by Eq. ( 1191) . 

Equations for the C-even and C-odd DD* matrix elements 

t± = t22 ± ^22) 

take the form: 

^^^ e) = -^pT 5 s(p-p)t 9 D ^_ ap _ pl) W) 

±9 ^)] d S D 3 (p-s-ap)D 2 (s) t ± 
which can be brought to the form of Eq. ( |26l) by the substitution 

tT(p, P ',E) = - nP l D i P) 5™5{p- P >) + ^ar(p,p'^) D2b,) 



A(p) v ^ ^ y A(p) * v *^' 7 A(p') ' 



